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ABSTRACT 

We present spectroscopic observations from the Spitzer Space Telescope of six carbon-rich 
AGB stars in the Sagittarius Dwarf Spheroidal Galaxy (Sgr dSph) and two foreground Galac- 
tic carbon stars. The band strengths of the observed C2H2 and SiC features are very sim- 
ilar to those observed in Galactic AGB stars. The metallicities are estimated from an em- 
pirical relation between the acetylene optical depth and the strength of the SiC feature. The 
metallicities are higher than those of the LMC, and close to Galactic values. While the high 
metallicity could imply an age of around 1 Gyr, for the dusty AGB stars, the pulsation peri- 
ods suggest ages in excess of 2 or 3 Gyr. We fit the spectra of the observed stars using the 
DUSTY radiative transfer model and determine their dust mass-loss rates to be in the range 
1.0-3.3 x 10 _8 Mq5t _1 . The two Galactic foreground carbon-rich AGB stars are located at 
the far side of the solar circle, beyond the Galactic Centre. One of these two stars show the 
strongest SiC feature in our present Local Group sample. 

Key words: circumstellar matter — stars: carbon — stars:AGB and post- AGB — stars: mass 
loss — galaxies: individual: Sagittarius dwarf Spheroidal — infrared: stars 



1 INTRODUCTION 

The Asymptotic Giant Branch (AGB) phase occurs during the late 
stages of the evolution of low- and intermediate-mass stars. This 
phase is characterised by intense mass loss and leads to the forma- 
tion of a circumstellar envelope made of gas and dust. The molec- 
ular composition of this envelope is dependent on the C/O abun- 
dance ratio. The CO molecule is very stable and unreactive, and 
if C/O > 1, all the oxygen is trapped in CO. The envelope (and 
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star) is then "carbon-rich"; amorphous carbon dominates the dust, 
while (after H2) CO and C2H2 dominate its gas. The"oxygen-rich" 
AGB stars are characterised by silicate dust and molecules such as 
SiO, OH, H2O. Third dredge-up brings the carbon produced by the 
triple-a reactions to the surface, increasing the C/O ratio and over 
time can change a star from oxygen-rich to carbon-rich. 

The mass loss from AGB stars is one of the main agents for 
the chemical evolution of galaxies. It expels the products of nuclear 
reactions in the core of the star into the interstellar medium. The 
mass loss from AGB stars contributes roughly half of all the gas 
recycled by stars (Maeder 1992), and up to nearly 90% of the dust 
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(Gehrz 1989). Mass loss from AGB stars is also one of the main 
sources of carbon in the universe, together with Wolf-Rayet stars 
and supernovae (Dray et al. 2003). 

The mechanisms driving this mass-loss process are not fully 
understood. It is thought to be a two-step process: pulsations from 
the star lead to the formation of dust, and then radiation pressure ac- 
celerates the dust grains to the escape velocity. The gas is then also 
expelled due to friction with the dust grains. The effect of metal - 
licity on the mass-loss rates from AGB stars has been discussed at 
some length. Some works proposed that mass-loss rates should be 
lower in metal-poor environments, as less dust is expected to form, 
so that radiation pressure would be less efficient in driving the mass 
loss (Bowen & Willson 1991; Zijlstra 2004). This hypothesis on 
the metallicity effect on mass loss has recently been tested using 
spectroscopy from Spitzer (e.g. Sloan et al. 2006, 2008; Zijlstra 
et al. 2006b; Lagadec et al. 2007; Matsuura et al. 2007; hereafter 
"SZLM" will refer to all five of these papers). Against early expec- 
tations, the dust mass-loss rates of carbon stars in metal-poor galax- 
ies from the Local Group appear to be similar to the ones measured 
in the Galaxy. In contrast, oxygen-rich stars in metal-poor environ- 
ments do appear to have lower dust mass-loss rates. This suggests 
that carbon is important in triggering the superwind (Lagadec & 
Zijlstra 2008). Theoretical models confirm that the mass-loss rates 
of carbon stars should not depend on metallicity (Wachter et al. 
2008). Mattson et al. (2008) argue that the pulsation energy of the 
star can drive a strong wind at low metallicity, but the dependence 
on chemistry suggests this is not the dominant effect within the ob- 
served metallicity range. 

The dusty circumstellar envelope surrounding AGB stars ab- 
sorbs the light from the central star and re-emits it in the in- 
frared. Furthermore, spectral signatures of dust and molecules ap- 
pear in this wavelength range. Spectroscopic observations of AGB 
stars in the thermal infrared are therefore vital for the study of 
the dusty envelopes around these stars. Atmospheric absorption 
prevents ground-based mid-infrared observations outside two win- 
dows around 10 and 20^tm. The Spitzer Space Telescope (Werner 
et al. 2004), with its high sensitivity and mid-infrared wavelength 
coverage, has proven to be a valuable tool for the study of dusty 
envelopes around extragalactic AGB stars. 

To study the mass loss from evolved stars at low metallicity, 
we are undertaking a Spitzer spectroscopic survey of mass-losing 
AGB stars in different Local Group galaxies. We have already 
presented results for the Magellanic Clouds and Fomax (SZLM; 
Matsuura et al. 2006; Groenewegen et al. 2007). Here we present 
Spitzer spectra of eight AGB stars in the direction of the Sagittarius 
Dwarf Spheroidal (Sgr dSph) galaxy^ This study aims at studying 
the circumstellar properties of these AGB stars. 



2 TARGET SELECTION 

We selected eight AGB candidates in the direction of the Sgr dSph 
based on their near-infrared colours. These stars were not spectro- 
scopically confirmed carbon stars before this study. The core of 
the Sgr dSph has a distance modulus estimated to be 17.02±0.19 
(Mateo et al. 1995) and is located behind the Galactic Bulge. It is 
being disrupted by the Galaxy so that its tidal tail surrounds the 
Milky Way. It contains several stellar populations. The dominant 

1 This galaxy is also known as the Sagittarius Dwarf Elliptical Galaxy 
(SagDEG), but it should not be confused with the Sagittarius Dwarf Irregu- 
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Figure 1. The Mk vs J — K diagram of the observed sample. Crosses 
represent the observed Sgr dSph targets and the large circle the foreground 
targets. Points represent the stellar population of Sgr dSph. Mk is the abso- 
lute magnitude and was calculated assuming a distance modulus estimated 
to be 17.02±0.19 (Mateo et al. 1995). 



one has a metallicity in the range [Fe/H]=— 0.4 to —0.7 and an 
age of 8. Oil. 5 Gyr (Bellazzini et al. 2006). A second population 
is younger and more metal-rich with [Fe/H]=— 0.25 for the most 
metal-rich objects (Zijlstra et al. 2006a). These populations span 
the range of metallicities observed in the Magellanic Clouds. 

The eight selected targets are a subsample of the stars pre- 
sented by Lagadec et al. (2008). TableQ] lists some characteristics 
of these stars. For convenience, we will use the names published 
previously. The stars were selected to span a wide range in J — K 
colour, i.e. a wide range of optical depth due to dust in their en- 
velopes. 

Fig. Q] displays the eight observed stars in an Mk vs J — if 
diagram. Asterisks represent the observed stars in the Sgr dSph. 
To show that the observed stars are all AGB stars, we overplot the 
distribution of Sgr field stars. Lagadec & Zijlstra (2008) explain 
the selection of these field stars. To verify that the observed stars 
belong to Sgr dSph, we performed radial velocity measurements: 
these indicate that six out of the eight observed stars belong to Sgr 
dSph but the remaining two are foreground stars (see Sec. |3.3l l. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Spitzer observations 

The observations were made with the InfraRed Spectrograph (IRS, 
Houck et al. 2004), on board the Spitzer Space Telescope. We used 
the Short-Low (SL) and Long-Low (LL) modules to cover the 
wavelength range 5-38^tm. The SL and LL modules are each di- 
vided in two spectral segments, together known as SL2, SL1 , LL2 
and LL1; a "bonus" order covering the overlap between the two 
modules is also available. The data reduction is similar to that de- 
scribed by Zijlstra et al. (2006b). The raw spectra were processed 
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through the Spitzer pipeline S 1 5 . We replaced the bad pixels by val- 
ues estimated from neighbouring pixels. The sky was subtracted by 
differencing images aperture by aperture in SL and nod by nod in 
LL. We used the software tools available in SPICE (the Spitzer IRS 
Custom Extractor) to extract the spectra. The flux calibration made 
use of the reference stars HR 6348 (KO III) in SL and HR 6348, 
HD 166780 (K4 III) and HD 173511 (K5 III) in LL. The spec- 
tra were individually extracted from the individual images. Both 
nods in both apertures were then joined simultaneously, recalculat- 
ing the errors in the process by comparing the nods. The different 
nods were averaged, using the differences to estimate the errors. 
The different spectral segments were combined using scalar mul- 
tiplication to eliminate the discontinuities due to flux lost because 
of pointing errors. The different segments were also trimmed to re- 
move dubious data at their edges. We also retained the data in the 
bonus order where it was valid. These steps resulted in a standard 
wavelength calibration accuracy of 0.06 /im in SL and 0.15 /im in 
LL. 

Fig-IHpresents the spectra of the observed AGB stars, ordered 
by their [6.4]-[9.3] colours (see Sec. 5). The molecular bands and 
dust features discussed below identify all eight objects as carbon- 
rich stars. 



3.2 Near-infrared photometry 

Observations in the near-infrared are important to obtain a reliable 
estimate of the luminosity. Due to the variability, these observations 
are best obtained close to the same epoch as the spectra. Multi- 
epoch JHKL photometry was obtained at the Australian National 
University (ANU) 2.3-m telescope at Siding Spring Observatory 
(SSO) in Australia. The filters used were centred at 1.28/im (J), 
1.68/im (H), 2.22/im (K) and 3.59/^m (L). Groenewegen et al. 
(2007) describe the observations. TableQ] presents the measured 
JHKL magnitudes at the epoch of the Spitzer observations. 

The multi-epoch observations, taken before the Spitzer ob- 
servations, were used to study the near-infrared variability of the 
sources. Light curves were obtained by fitting a sine wave to the 
J\~-band data. Fig.[3]displays these light curves. 



3.3 Radial velocities 

Radial velocities were determined for six of the eight carbon stars 
using optical spectra obtained with the Dual Beam Spectrograph 
on the 2.3-m SSO telescope. The spectra have a resolution of 0.48 
A/pixel for four objects and 3.7 A/pixel for the remaining two 
objects (see TableQ. The higher resolution spectra were cross- 
correlated with the local carbon star X Vel for which a spectrum 
at a resolution of 0.5 A/pixel was also obtained. A radial veloc- 
ity of —5.4 kms" 1 was adopted for X Vel (determined from high- 
resolution echelle spectra which were cross correlated against the 
radial velocity standard a Cet with a heliocentric radial velocity of 
—25.8 km s _1 ). Cross-correlation of the lower-resolution spectra of 
the remaining two stars was used to obtain their radial velocities. 
Table|2]gives the final heliocentric radial velocities and their errors 
(as given by the IRAF task FXCOR). 

Fig. ID presents a histogram of the radial velocities of carbon 
stars in the direction of Sgr dSph from Ibata et al. (1997). There 
are obviously two groups of stars: those with heliocentric radial 
velocities of 100 < vi lc ii (km/s) < 190 which belong to the Sgr 
dSph; and those with Vh c ii (km/s) < 50 which belong to the Milky 
Way Galaxy. Two stars in our sample, Sgr 02 and Sgr 22, clearly 
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Figure 4. Top panel: Radial velocity distribution of stars in the direction of 
the Sgr dSph galaxy (from Ibata et al. 1997). Bottom panel: radial velocity 
distribution of six of our targets. Two of these stars, Sgr 02 and Sgr 22, are 
clearly foreground members of the Galaxy. 



Table 2. Heliocentric radial velocities for six of the observed targets as 
measured from observations at SSO. 



Target 



Resolution 
A/pixel 



wavelength range 

°A 



Vhelio 

km/s 



km/s 



Sgr 02 


0.48 


8100-8900 


-22.8 


1.1 


Sgr 03 


3.7 


6800-9200 


128.1 


11.1 


Sgr 09 


0.48 


8100-8900 


144.9 


1.6 


Sgr 15 


3.7 


6800-9200 


136.4 


12.2 


Sgr 22 


0.48 


8100-8900 


-73.1 


2.1 


Sgr 29 


0.48 


8100-8900 


106.3 


2.5 



belong to the Milky Way Galaxy, while the other stars are members 
of the Sgr dSph. The properties of the two foreground stars are 
discussed in Sec. 19. 31 

We used other distance estimates to ascertain the Sgr dSph 
membership of the two stars without radial velocity measurements 
(see SeclTJl. 



4 DESCRIPTION OF THE SPECTRA 

The IRS spectra of the eight observed stars (Fig. [2j show both 
dust and molecular emission. All the observed stars are carbon- 
rich and have spectra typical of AGB stars. The spectra of all stars 
clearly show absorption features from C2H2 at 7.5 and 13.7 
(Fig. HJ. SiC dust produces the 11.3 /im feature detected in all the 
spectra. Some of them (Sgr02, Sgr03, Sgr07 and Sgr 18) show a 
broad emission feature around 30 fim attributed to MgS (Hony et 
al., 2002). This feature might also be present in Sgr 09 and Sgr 22. 
A drop is observed for all stars at the blue edge of the spectra. It is 
due to several molecules, most notably C3 and CO (Zijlstra et al. 
2006b). 

Sgr 22 displays a very strong SiC emission feature, stronger 
than any carbon-rich AGB star of which we are aware, including 
the Galactic sample observed by the Infrared Space Observatory 
and described by Sloan et al. (2006) or the Local Group galaxies 
observed by Spitzer (see upper panel of Fig.[8j. 

A weak feature attributed to the stretching vibration of a car- 
bonyl group (X-CO) has been observed in some LMC AGB stars 
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ure 2. Spitzer/IRS spectra of the eight observed carbon stars, ordered by dust temperature. The main dust and gas features are labelled on the spectrum of 
02. The dashed lines represent the SED fits obtained with DUSTY (Sec. [8) 
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Figure 3. Light curves of the eight observed carbon stars. These were obtained using a sinusoidal fit to the ii"-band magnitudes. The vertical lines show the 
dates of our Spitzer observations. 
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Table 1. Observed Sgr dSph targets: adopted names, coordinates, photometry and distance. J, H, K and L are taken from near-simultaneous measurements 
at SSO. The periods were obtained by sine-fit to the A"-band observations. The distance was estimated using the method described in Sec. [7] 



Adopted name 


2MASS name 


IRAS name 


RA 


Dec 


J 


H 


K 


L 


P 


phase 


D 












(J2000) 


mag 


mag 


mag 


mag 


d 




kpc 


Sgr 02 


18414350 


-3307166 


18384- 


-3310 


18 41 43.50 


-33 07 16.6 


11.095 


9.421 


8.043 


6.501 


301 


0.65 


17.1 


Sgr 03 


18443095 


-3037098 


18413- 


-3040 


18 44 30.96 


-30 37 09.8 


13.360 


11.436 


9.879 


8.047 


446 


0.36 


31.6 


Sgr 07 


18465160 


-2845489 


18436- 


-2849 


18 46 51.60 


-28 45 48.9 


15.464 


13.121 


10.944 


8.320 


512 


0.06 


34.1 


Sgr 09 


18514105 


-3003377 






18 51 41.05 


-30 03 37.7 


11.753 


10.263 


9.124 


7.913 


370 


0.40 


30.3 


Sgr 15 


18584385 


-2956551 


18555- 


-3001 


18 58 43.85 


-29 56 55.1 


12.862 


10.921 


9.394 


7.496 


417 


1.00 


28.0 


Sgr 18 


19043562 


-3112564 


19013- 


-3117 


19 04 35.62 


-31 12 56.4 


15.535 


12.753 


10.734 


8.171 


485 


0.42 


35.1 


Sgr 22 


19103987 


-3228373 


19074- 


-3233 


19 10 39.87 


-32 28 37.3 


11.315 


9.454 


7.959 


6.196 


370 


0.02 


16.0 


Sgr 29 


19485065 


-3058317 






19 48 50.65 


-30 58 31.9 


12.909 


11.014 


9.522 


7.838 


339 


1.00 


31.5 



(Zijlstra et al. 2006b) at 5.8/im. This feature is present in the spec- 
tra of Sgr 09 and Sgr 18, and possibly also in Sgr 02, Sgr03, Sgr 07 
and Sgr 15. 



5 COLOURS AND BAND STRENGTHS 

To determine colours of the observed stars, we applied the so-called 
"Manchester System" (Sloan et al. 2006; Zijlstra et al. 2006b), us- 
ing four narrow bands selected to represent the continuum at differ- 
ent wavelengths. Using this method, we derive two colours. Zijlstra 
et al. (2006b) showed that the [6.4] — [9.3] colour is a good esti- 
mate of the dust optical depth, while the [16.5] — [21.5] colour pro- 
vides an estimate of the dust temperature. The [6.4] — [9.3] colour 
shows a linear relation with the measured mass-loss rates (Groe- 
newegen et al. 2007; Matsuura et al. 2007; Sloan et al. 2008). Table 
[3] lists the blackbody temperature derived from this [16.5] — [21.5] 
colour. To test that our choice of continuum colours is reasonable, 
we plotted [16.5] — [21.5] vs. [6.4]-[9.3] (Fig.©. For five of the 
observed stars, the continua are consistent. But three stars (Sgr 02, 
Sgr 03 and Sgr 22) are outliers. Sgr 02, Sgr 03 are particularly red 
at[16.5]-[21.5]. 

For all the observed stars, we measured the strength of the 
main features (Fig. [2}. The continuum underlying each feature is 
defined using small wavelength ranges on the blue and red sides 
of the feature. The continuum is then defined using a straight 
line between the blue and red continuum values. Because the red 
edge of the MgS feature (Sec. @ is outside of the IRS spectral 
range, we used a blackbody with the temperature derived from the 
[16.5] — [21.5] to extrapolate the continuum under this dust feature. 
After the definition of the continuum we determined the strengths 
of the features in two different ways. For the dust features, we sim- 
ply measured the ratio between the integrated feature flux and the 
continuum (F/C). For the gas features, seen in absorption, we mea- 
sured the equivalent widths. We also measured the central wave- 
length of the SiC feature defined as the wavelength at which the 
flux on the blue side of the feature equals the flux on the red side. 
Table[3]lists the measured strengths and central wavelengths of the 
observed features. 



6 CIRCUMSTELLAR PROPERTIES 
6.1 Gas 

As mentioned in Sec. [4] the main gas absorption features observed 
in our spectra at 7.5/im and 13.7/im are due to C2H2. We have 
shown (SZLM, Matsuura et al. 2006) that in metal-poor environ- 
ments, the C2H2 absorption becomes stronger. This can be ex- 
plained by the fact that at low metallicity, the initial oxygen abun- 
dance is smaller than for Galactic stars, and the dredge-up of carbon 
during the AGB leads to higher C/O ratios. 

Fig.|7]shows the equivalent width of these features as a func- 
tion of the [6.4] — [9.3] colour for the current sample in comparison 
to the stars from other Local Group galaxies. The plot relates the 
strength of the molecular band to the optical depth (or dust mass- 
loss rates) of the envelopes of the observed stars. The C2H2 equiv- 
alent widths of the Sgr dSph sample are in the lower range of the 
observed strengths. They are similar to or a little higher than the 
values found in Galactic stars of similar [6.4] — [9.3], but are gener- 
ally weaker than found in the SMC, LMC and Fornax. 



6.2 Dust 

As mentioned in Sec. [4] the main dust features observed in the 
spectra are due to SiC and MgS. The featureless continuum arises 
primarily from emission from amorphous carbon emission. Fig. [8] 
shows the strength of the SiC and MgS features as a function of the 
[6.4] -[9.3] colour. 

The properties of the SiC feature have been extensively stud- 
ied previously (e.g. Speck et al. 2005, Leisenring et al. 2008). La- 
gadec et al. (2007) have shown that the relationship between the 
strength of the SiC feature and the optical depth of the envelope 
varies according to the metallicity of the host galaxy, using a sam- 
ple of stars in the SMC, LMC and the Galaxy. Fig.[8]illustrates the 
trend of increasing SiC feature strength for the Galactic stars when 
[6.4] — [9.3] Ss0.5 and a gradual decrease for redder stars. Such a 
trend is observed for carbon stars in the LMC, but the inflexion 
occurs for redder colours, around ~1. Very few SMC and Fornax 
stars with [6.4] — [9.3] > 1 have been observed, but such stars in 
metal-poor galaxies appear to follow the same trend as LMC stars 
and have much weaker SiC features at a given optical depth. One 
SMC star, GM 780, is an outlier with a very strong SiC emission 
feature. The stars from the Sgr dSph are located at similar positions 
as the Galactic ones, but with a SiC feature slightly weaker than 
those of Galactic stars with the same [6.4]-[9.3] colour. As Si is not 
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Table 3. Colours measured using four narrow carbon stars continuum bands, strength of the molecular and dust features, in terms of either the equivalent 
width in microns, or the integrated fiux-to-continuum ratio, F/C (Sect. [5)- The last column gives the continuum (black-body) temperature, derived from the 
[16.5] — [21.5] colour listed in TableQ] 



target [6.4]-[9.3] [16.5]— [21.5] EW (7.5 pm) EW (13.7 /mi) F/C(SiC) A c (SiC) F/C (MgS) T( 

(mag) (mag) (pm) (pm) (pm) (K) 



Sgr 02 


0.375 ±0.012 


0.270 ±0.010 


0.148± 0.002 


0.035 ± 0.005 


0.267± 0.005 


11.27 ±0.03 


0.639 ± 0.012 


515± 16. 


Sgr 03 


0.555 ±0.004 


0.292 ±0.009 


0.1 15± 0.005 


0.020 ± 0.008 


0.226± 0.007 


11.07 ±0.05 


0.401 ±0.012 


480.± 12. 


Sgr 07 


0.757 ±0.004 


0.187 ±0.015 


0.078± 0.005 


0.036 ± 0.002 


0.157± 0.004 


11.17 ±0.04 


0.180 ±0.018 


716.±48. 


Sgr 09 


0.319 ±0.016 


0.093 ±0.018 


0.237± 0.009 


0.064 ±0.011 


0.278± 0.008 


11.14 ±0.04 


0.050 ± 0.026 


1385.± 225. 


Sgr 15 


0.532 ±0.006 


0.152 ±0.014 


0.061 ±0.004 


0.018 ± 0.005 


0.222± 0.005 


11.11 ±0.03 


0.172 ±0.018 


864.± 66. 


Sgr 18 


0.821 ±0.002 


0.188 ±0.013 


0.068± 0.004 


0.034 ± 0.003 


0.158± 0.003 


11.18 ±0.03 


0.234 ±0.015 


710.± 41. 


Sgr 22 


0.636 ±0.005 


0.089 ±0.016 


0.093± 0.002 


0.023 ± 0.002 


0.410± 0.006 


11.20 ±0.02 


0.055 ± 0.020 


1457.± 225. 


Sgr 29 


0.362 ±0.010 


0.105 ±0.017 


0. 13 1± 0.004 


0.044 ± 0.010 


0.1 89± 0.008 


11.22 ±0.06 


0.001 ± 0.021 


1230.± 166. 



produced in AGB stars, this slight difference appears to indicate 
that the Si abundance in the observed Sgr dSph stars is between the 
Si abundance of Fornax/SMC/LMC and the Galaxy. 

Fig. [9] shows the strength of the MgS feature as a function of 
the dust temperature. As already observed for AGB stars in other 
Local Group galaxies (Sloan et al. 2006, Zijlstra et al. 2006b, La- 
gadec et al. 2007), this feature tends to appear mostly in the en- 
velopes of stars with cool dust. In general, the formation process of 
MgS begins around 600 K and is complete around 300 K. 



7 DISTANCE ESTIMATES 

To check whether the observed stars are members of Sgr dSph, we 
estimated their distances using two methods. The distance to Sgr 
dSph being well established, our sample allows us to test different 
distance estimation methods. 

The first one uses the infrared colours of the observed stars, 
following the method described by Sloan et al. (2008). They have 
shown, from a sample of carbon stars in the Magellanic Clouds, 
that: 

M K = -9.18 + 0.395(J - #) (1) 
and 

A /«. •,..-. = a »([6-4] - [9.3])* (2) 

where ao = —8.81, ai = —6.56 and a,2 = 2.74. Using these 
relations and the infrared colours of the observed carbon stars, we 
can obtain two estimates of their distances. 

The second method makes use of the period-luminosity rela- 
tion for carbon Miras as described by Feast et al. (2006): 

M bo i = -2.54 x logP + 2.06 (3) 

We converted all of the photometry to the SAAO system and 
dereddened the near-infrared data as described by Lagadec et al. 
(2008). The reddening was estimated using the extinction maps 
by Schlegel et al. (1998). For the observed stars, Aj~0.1 and 
Ax~0.04. Using standard conversion factors, we estimated that 
A6/xm<0.01 and that the reddening was lower than a percent at 
longer wavelengths. We thus did not deredden our spectra. We then 
derived the bolometric magnitudes using the equation for bolomet- 
ric correction derived by Whitelock et al. (2006): 

BC K = ±0.972 + 2.9292 x (J - K) - 1.1144 x (J - K) 2 



Table 4. Distances estimated from the infrared colours and the period- 
luminosity relation. 



target 


D/h (kpc) 


DpL(kpc) 


Sgr 02 


17.1 


12.9 


Sgr 03 


31.6 


33.2 


Sgr 07 


34.1 


41.4 


Sgr 09 


30.3 


25.2 


Sgr 15 


28.0 


25.8 


Sgr 18 


35.1 


32.7 


Sgr 22 


16.0 


12.9 


Sgr 29 


31.5 


25.1 



+0.1595 x (J - K) 3 - 9.5689 103( J - K) 4 (4) 

Table [4] displays the mean values of the estimated distance 
for each star using the dereddened infrared colours and the period- 
luminosity relation. This confirms that the stars Sgr 07 and Sgr 18 
belong to Sgr dSph. 

Mateo et al. (1995) derived a distance for the central region 
of Sgr dSph of 25.35±2.06 kpc. A comparison of our two distance 
estimates indicates that the best agreement with the distance to Sgr 
dSph is obtained with the period-luminosity relation. Using the in- 
frared colours gives distances to the observed stars significantly 
higher than the well-determined distance to the galaxy. The dis- 
tance estimated using the period-luminosity relation gives results 
in perfect agreement for Sgr 09, Sgr 15 and Sgr 29. For the three 
other stars, this method gives larger distances. Variability will have 
some effect on the derived distances, but given that two of the three 
stars with large apparent distances are significantly redder in J — K 
than any of the others it would seem quite possible that they are un- 
dergoing obscuration events of the type described by Whitelock et 
al. (2006) for Galactic Miras and which are also seen in at least one 
of the AGB stars in Fornax (Whitelock et al. 2009). Photometric 
monitoring over long time scales of these stars would be needed to 
confirm that hypothesis, as stars with such behaviours show varia- 
tions over periods of a few years. 



8 DUST MASS-LOSS RATES 

Lagadec et al. (2008) measured the dust mass-loss rates of the stars 
presented in this work. These dust mass-loss rates (M co i), mea- 
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Table 5. Dust mass-loss rates for the observed stars and parameters obtained from our DUSTY models. M col and M dusty are the dust mass-loss rates (in 
MQyr - 1 ) from Lagadec et al. 2008 and our DUSTY models, respectively. 



Target Luminosity A^DUSTY ^eff T; n SiC/AMC r (0.55/xm) M co \ V cxp 

(L ) (10- 8 M Q yr- 1 ) (K) (K) (lO^Meyr" 1 ) km/s 



Sgr02 


11469 


2.59 


2800 


1200 


0.1 


4.67 


0.86 


28 


Sgr03 


4529 


1.29 


2800 


1200 


0.1 


9.56 


0.84 


18 


Sgr07 


4483 


2.52 


2800 


1000 


0.08 


9.64 


0.97 


14 


Sgr09 


6652 


1.02 


2800 


1200 


0.1 


2.22 


0.27 


27 


Sgrl5 


8632 


2.73 


2800 


1200 


0.1 


7.00 


0.93 


23 


Sgr 18 


5056 


3.39 


2800 


1000 


0.08 


13.1 


2.60 


13 


Sgr22 


13321 


3.78 


2800 


1200 


0.2 


7.11 


1.30 


26 


Sgr29 


9152 


2.04 


2800 


1500 


0.1 


5.89 


0.96 


35 



sured from near and mid-infrared colours are presented in Table 

HI 

The Spitzer spectra, together with simultaneous near-infrared 
photometric measurements give us an opportunity to better model 
the dust emission from these stars. We modeled the spectral energy 
distributions (SEDs), defined by the JHKL flux and the Spitzer 
spectra, using the radiative transfer code DUSTY (Ivezic & Elitzur 
1997). Our primary objective is to fit the dust continuum in order to 
estimate the optical depths and the mass-loss rates. DUSTY solves 
the 1-D problem of radiation transport in a dusty environment. For 
all our models, we assume that the irradiation comes from a point 
source (the central star) at the centre of s spherical dusty envelope. 
The circumstellar envelope is filled with material from a radiatively 
driven wind. All the stars are carbon-rich and display SiC in emis- 
sion. We thus assumed that the dust composition of the envelope 
is a mixture of amorphous carbon and SiC. Optical properties for 
these dust grains are taken the work of Hanner (1988) and Pegourie 
(1988) for amorphous carbon and SiC, respectively. The grain-size 
distribution is assumed to be a typical MRN distribution, with a 
grain size a varying from 0.0005 to 0.25/im distributed according 
to a power law with n(a)oca~ q with q=3.5 (Mathis et al. 1977). The 
outer radius of the dust shell was set to 10 3 times the inner radius; 
this parameter has a negligible effect on our models. 

To model the emission from the central star, we used a hy- 
drostatic model including molecular opacities (Loidl et al. 2001, 
Groenewegen et al. 2007). We used those hydrostatic models as an 
illuminating source only, we did not attempt to model the molecu- 
lar absorption features. The C/O ratio is assumed to be 1.1 for the 
hydrostatic models, which is a typical value for Galactic AGB stars. 
This ratio is not well known in the observed stars: the C/O ratio of 
AGB stars appears to be higher in metal-poor environments than in 
the Galaxy (Matsuura et al. 2002, 2005). But varying the C/O ratio 
mainly affects the depth of the molecular features (such as CO and 
C2H2). The radiative transfer model fits the dust emission, and the 
C/O ratio has little effect on the model. We also performed some 
model calculations assuming that the central star emits as a black- 
body, but these did not provide satisfactory fit to the near-infrared 
data, due to the absence of molecular absorption. The free parame- 
ters of the models are the dust temperature at the inner radius, the 
mass ratio of SiC to amorphous carbon dust and the central star 
effective temperature, T e fj. As one output, DUSTY provides an es- 
timate of the terminal outflow velocity and the mass-loss rate for a 
luminosity of 10 4 Lq. As the luminosity can be determined by scal- 
ing the emerging spectrum, we can then recalibrate the mass-loss 



rates accordingly, assuming a distance of 25 kpc for the six stars in 
Sgr dSph, and 17.1 and 16.0 kpc for Sgr02 and Sgr22 respectively. 

DUSTY gives total (gas+dust) mass-loss rates assuming a gas- 
to-dust ratio of 200, but the dust is the dominant constraint on the 
fitting to the SED. Therefore, we have divided the resulting mass- 
loss rates by the gas-to-dust ratio and report only the dust mass- 
loss rate in Table [5] This assumes that the gas and dust expansion 
velocities are the same. Note that we had to make a number of 
assumptions in our models. For example, we assume that the enve- 
lope and the dust grains are spherical. This cannot be the case in 
reality, so our fits cannot be perfect. However, the models achieve 
our primary aim, estimating the dust mass-loss rates. These values 
are of the same order of magnitude as those derived by Lagadec et 
al. (2008). The values found by Lagadec et al. (2008) are generally 
a factor of two or three smaller than the present ones. 

For all the stars, the best fit is obtained with T c ff=2800K. The 
best model for most of the stars was obtained using a temperature at 
the inner radius of 1200 K and a mixture of 10% of SiC and 90% of 
amorphous carbon. For two stars, we needed Ti n = 1000 K to model 
the red part of the SED, and for another one we used Ti n = 1500 K. 
The very strong SiC feature observed in the spectrum of Sgr 22 is 
very well fitted using 20% of SiC. Sgr 07 and Sgr 18 have weaker 
SiC features, which can be fitted using a SiC mass fraction of 8%. 

We confirm the general correlation found by Groenewegen 
et al. (2007) and Matsuura et al. (2007) between the mass-loss 
rates and the [6.4] — [9.3] (FigllOt- But the stars we observed in 
the Galaxy and in Sgr dSph are on average offset towards higher 
mass-loss rates for the same colour, by about a factor of two. The 
two stars in Fomax show a (smaller) offset in the same direction. 
All of these mass-loss rates of these stars were estimated from our 
DUSTY models, while the ones in the Magellanic Cloud stars were 
estimated using another radiative transfer model (Groenewegen et 
al., 2007). The models by Groenewegen et al. (2007) were made as- 
suming a constant expansion velocity of 10 km s _1 while DUSTY 
calculates an expansion velocity for each model. These calculated 
velocities are generally a factor of 2-3 higher than the ones assumed 
by Groenewegen et al., suggesting that the difference in mass-loss 
rates arise from different assumptions made by the models, rather 
than an intrinsic differences between the stars. The dust mass-loss 
rates determined by Lagadec et al. (2008) for the Sgr dSph stars 
nicely follow the relation with [6.4] — [9.3] colour found in previ- 
ous work. 
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Figure 5. Spitzer spectra around 13.7/xm of the observed stars, showing 
the C2H2 absorption. The two grey bands represent the continuum used to 
measure the band strengths. 



9 DISCUSSION 

9.1 Estimating metallicities 

The survey of carbon-rich AGB stars in the Local Group has 
yielded Spitzer IRS spectra of a large sample of stars in galaxies 
covering a range of metallicities. 

The metallicity of extragalactic AGB stars is normally taken 
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Figure 6. [16.5] — [21.5] vs. [6.4] — [9.3] colour diagram. Asterisks and open 
circles represent the stars we observed. The solid line and filled squares 
represent a series of DUSTY models, as described by Zijlstra et al. (2006a). 
The optical depth at 1 fj,m is indicated for each model 

as that of the underlying dominant stellar population of the galaxy. 
Individual metallicities of carbon stars are especially difficult to 
determine. Generally, stellar metallicities are determined from the 
strengths of fine-structure lines from heavy elements in optical or 
near-infrared spectra. The continua of AGB stars are very uncer- 
tain, due to the presence of multiple, overlapping molecular bands 
in their envelopes, e.g. CO and CN. To model this continuum, a 
hydrostatic stellar atmosphere model is required, with a known sur- 
face gravity, C/O abundance ratio, and effective temperature. Such 
models have been made to determine the metallicity of Local Group 
AGB stars, but are very time-consuming, and have been applied to 
only a small number of stars (de Laverny et al. 2006). 

The SiC feature is expected to depend on metallicity, through 
the Si abundance. It may be possible to use the strength of this fea- 
ture in our Spitzer spectra to estimate the metallicity of the stars 
in our sample. Simultaneously, the C2H2 abundance depends on 
the amount of free carbon, [(C— 0)/0] (Lagadec & Zijlstra 2008), 
which is expected to increase (on average) with decreasing metal- 
licity (e.g. Matsuura et al. 2005). Thus the Spitzer spectra provide 
two means of determining metallicity. 

Fig.ll llshows the equivalent width of the C2H2 7.5 /im feature 
versus the SiC flux-over-continuum ratio, for all the Local Group 
carbon stars we observed. A clear separation is seen between the 
stars from different galaxies, i.e. different metallicities. The Galac- 
tic stars, which are expected to be the most metal-rich, have the 
largest SiC/C ratio for a given C2H2 7.5 equivalent width. The 
LMC stars are spread out between the SMC and Galactic stars. The 
diagram shows a good (average) separation between the galaxies, 
in order of expected metallicity. Stars from each individual galaxy 
show a significant spread, which may be due to a spread in metal- 
licity within a galaxy, or due to an evolutionary spread, as the stars 
will show an increasing carbon dredge-up over time. But for the full 
carbon-star population, this diagram appears to provide a metallic- 
ity indicator. 

We can apply this diagram to the stars observed in this paper. 
The two galactic foreground stars fall within the general population 
of Galactic stars, with one star (Sgr 22) showing an unusually high 
SiC/C ratio which can be interpreted as arising from a high metal- 



10 E. Lagadec et al. 



0.15 



3 
a. 



6 
a. 



o.i 



O Galaxy 
■ SMC 
• LMC 



0.10 - S> 



0.05 



n o°e 



• • • 



°4> ■. * 



° ■ A 



0.2 



O o . 



0.4 



0.3 



A Sgr dSph 
O Galaxy 



• LMC 

A Fornax 



* 0.2 



4& 

o 
o o 

A 



u a • 



0.1 :* d 



0.0 



£ • ■ 



■ .i ° 

■ ■ ■ ■ %m* • ■ 



■ a 

% ■ ■ 



0.0 0.5 1.0 1.5 

[6.4]-[9.3] 

Figure 8. The strength of the SiC features as a function of the [6.4]— [9.3] 
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(Sec.[5j. Comparison data are from the literature are as in Fig[7] 



0.0 L. 
0.0 



0.5 1.0 1.5 

[6.4]-[9.3] 



2.0 



Figure 7. The equivalent width of the C2H2 features at 13.7^m (upper 
panel) and 7.5/im (lower panel) as a function of [6.4] — [9.3] colour. The 
plot combines stars from the present sample and stars from previous Lo- 
cal Group observations: open triangles represent the current Sgr dSph sam- 
ple, filled squares the SMC stars from Sloan et al. (2006) and Lagadec et 
al. (2007) samples, open circles the ISO Galactic sample defined by Sloan 
et al. (2006), filled circles the LMC sample of Zijlstra et al. (2006b) and 
Leisenring et al. (2008), and triangles stars from the Matsuura et al. (2007) 
Fornax sample. 



licity. Surprisingly, the stars in the Sgr dSph are found between the 
Galactic and LMC stars, with an indicated metallicity higher than 
that of the LMC, although still sub-solar. Based on this diagram, the 
observed Sgr dSph carbon stars are the most metal-rich population 
after the Galactic stars. 
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Figure 9. Strength of the MgS feature as a function of dust temperature 
estimated from the [16.5] — [21.5] colour. Literature data as in Fig|7] 



9.2 Metal-rich carbon stars in the Sgr dSph 

The dominant giant branch population in the Sgr dSph is believed to 
have a metallicity ([Fe/H]) ~-0.55 or less (Whitelock et al. 1996, 
Dudziak et al. 2000, Carraro et al. 2007) which is lower than that 
of the LMC. The observed high metallicity implied for the dusty 
carbon stars in our Sgr dSph sample is therefore unexpected. 

Evidence for a more metal-rich stellar population in the Sgr 
dSph was first detected by Bonifacio et al. (2004), who, based on 
spectroscopy of ten giants, found the dominant population to have 
[Fe/H] = —0.25, much higher than what is found from typical 
AGB stars. Zijlstra et al. (2006a) and Kniazev et al. (2008) show 
that one of the four planetary nebulae in the galaxy has the same 
high metallicity, confirming the existence of this population even if 
it does not dominate. Chou et al. (2007) find a median metallicity 



in the core of the Sgr dSph of [Fe/H] = -0.4, similar to the LMC 
and again a little higher than typical AGB stars. 

The dusty carbon stars studied in this paper appear to be com- 
paratively metal-rich. It is therefore natural to suggest that they 
could be an AGB counterpart of the population whose age has been 
estimated at 1 Gyr (Bonifacio et al. 2004), the formation of which 
may have been triggered by the passage of the dwarf galaxy through 
the Galactic disc. However, their pulsation periods, which range 
from 310 to 512 days, imply considerably older ages. For example 
there are three carbon Miras which have periods around 490 days 
(Nishida et al. 2000) in Magellanic Cloud Clusters which have ages 
about 1.6 Gyr (Mucciarelli et al. 2007a,b; Glatt et al. 2008), whilst 
van Loon et al. (2003) suggest that the cluster KMHK 1603, which 
contains a carbon-rich Mira with a period of 680 days, has an age of 
0.9-1.0 Gyr. The kinematics of Galactic carbon-rich Miras (Feast et 
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Figure 10. Dust mass-loss rates as a function of the [6.4] — [9.3] colour. The 
comparison data are described in Fig. [7] 



Figure 11. Strength of the C2H2 7.5/xm feature as function of the strength 
of the SiC feature. The comparison data are as described in Fig[7j 



al. 2006) is also consistent with increasing periods going with de- 
creasing age. The periods of the Sgr dSph Miras would therefore 
imply ages of 2 to 3 Gyr, or perhaps a little greater. More work is 
required to reconcile these observations. 

One can wonder why a relatively small galaxy would show 
the most metal-rich stellar population in the Local Group after the 
two main spirals (e.g. Zijlstra et al. 2006a). In general the dwarf 
spheroidals show a much steeper age-metallicity relation than do 
dwarf irregulars of the same luminosity. Dwarf spheroidals are 
found around large galaxies while the dwarf irregulars are more 
isolated systems, showing the effect of environment on enrichment 
history. It is also of interest that of the satellites of the Milky Way, 
only the Sgr dSph has reached such high metallicity. Fornax, other- 
wise a very similar system, contains carbon stars with more metal- 
poor characteristics (Matsuura et al. 2007). As the Sgr dSph is the 
nearest of the satellites (Fornax is around 130 kpc away), the effects 
of the Galaxy seem important. 

There are two possible explanations for the very high metallic- 
ities. First, the stripping of the interstellar gas by a passage through 
the Milky Way may have allowed fast chemical evolution through 
stellar mass loss and supemovae. Second, the system may have cap- 
tured gas from the Milky Way. The second option may have some 
support from the fact that the high metallicity is the same as that of 
the Galactic disk at the orbital distance of the Sgr dSph. 



9.3 Two foreground carbon stars 



As discussed in Sec. 13. 31 two stars in our sample, Sgr 02 and Sgr 22, 
belong to the Milky Way Galaxy while the other stars are members 
of the Sgr dSph. These stars have been selected through their J 
and K colours and were not spectroscopically confirmed carbon 
stars before these observations. The Spitzer spectra clearly show 
that these stars are carbon-rich AGB stars. Their distance is esti- 
mated to be 17.1 kpc and 16.0 kpc for Sgr 02 and Sgr 22, respec- 
tively, from their infrared colours (see Sec. [7}. Sgr 22 displays the 
strongest SiC feature among all the stars from our present Local 
Group sample (including the Galaxy). The other foreground star, 
Sgr 02, has a weaker SiC feature than Sgr 22, but Fig. [8] indicates 



that its SiC feature strength is similar to the ones observed in the 
Galaxy rather than in more metal-poor galaxies. 

The two Galactic stars are on line of sight to the Sgr dSph, 
and they are located on the far side of the Galaxy, close to the solar 
circle. There are very few carbon stars known in the inner Galaxy: 
carbon star distribution begins around the solar circle and extends 
outwards (Le Bertre et al. 2001). The two stars are therefore located 
around the inner edge of the carbon star distribution and should be 
similar to carbon stars in the Solar neighbourhood. However, their 
height above the Galactic plane is 3.6 and 4.9 kpc respectively. This 
is well beyond what is expected from the disc stars. There is no 
clear explanation, and it would be very interesting to establish the 
metallicity of other halo carbon stars, some of which are known to 
be Mira variables (e.g. Totten et al. 2000; Mauron 2008). The high 
metallicity suggests the stars escaped from the disc. Whether the 
escape process was triggered by the Sgr dSph passage is not clear. 



10 CONCLUSIONS 

We have presented a Spitzer spectroscopic survey of six carbon-rich 
AGB stars in Sgr dSph and two serendipitously discovered Galac- 
tic carbon stars. We used the low-resolution short-low and long-low 
mode of the IRS spectrometer to obtain mid-infrared spectra cov- 
ering the range 5-38/im. These observations enabled us to study 
dust and gas features from the circumstellar envelopes around these 
stars. We observed molecular absorption bands due to C2H2 at 7.5 
and 13.7/iin. A weak absorption feature at 14.3/im, due to HCN, is 
observed in two stars. A weak absorption feature at 5.8/im due to 
a carbonyl group (X-CO) is observed around two stars and might 
be present in the spectra of three others. The continuum of all the 
observed stars is due to emission from amorphous carbon, which 
does not display any spectral features. All the stars show a SiC fea- 
ture at 11.3jum. A broad dust emission emission due to MgS is also 
observed for six of the stars around 30/im. 

Radial velocity measurements show that two of our sample 
are actually carbon stars. The remaining six stars in our sample 
are members, as confirmed with radial velocities, infrared colour- 
magnitude relations, and/or period-luminosity relations. One of the 
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two Galactic carbon stars displays the strongest SiC feature ever 
observed. Both Galactic stars are certainly carbon stars located at 
the far side of the solar circle, but they have unexpectedly large 
distances from the Galactic plane. 

We fitted the SEDs of all the observed spectra using the radia- 
tive transfer code DUSTY. The estimated dust mass-loss rates are 
found to be in the range 1.0-3.3 x 10~ 8 M Q yr~\ 

The observed strengths of C2H2 and SiC are very similar to 
the ones observed for Galactic carbon stars. Stronger C2H2 and 
weaker SiC feature were expected in this metal-poor dwarf galaxy. 
We have shown that the strength of these features depends on the 
metallicity and that the observed stars have metallicities close to 
Galactic values. This result is unexpected, as the usual measure- 
ments of the metallicity ([Fe/H]) of the Sgr dSph range from —0.4 
to —0.7. The enhanced metallicity of the observed carbon stars 
indicates that the interstellar medium in the Sgr dSph has been 
strongly enriched. 
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